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Nonlinear current response of one- and two-band superconductors 
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We have calculated the nonlinear current of a number of single band s-wave electron-phonon 
superconductors. Among issues considered were those of dimensionality, strong electron-phonon 
coupling, impurities, and comparison with BCS. For the case of two bands, particular attention is 
paid to the role of anisotropy, the integration effects of the off-diagonal electron-phonon interaction, 
as well as inter- and intraband impurities. For the specific case of MgB2, we present results based 
on the known microscopic parameters of band theory. 
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I. INTRODUCTION 



With the advent of high temperature superconductiv- 
ity in the cuprates and the possibility of exotic gap sym- 
metry including nodal behavior, a renewed effort to find 
novel experimental probes of order parameter symme- 
try has ensued. One result of this effort was the pro- 
posal by Sauls and co-workers0 to examine the non- 
linear current response of d- wave superconductors. They 
showed that a nonanalyticity in the current- velocity rela- 
tion at temperature T — is introduced by the presence 
of nodes in the order parameter. One prediction was 
that an anisotropy should exist in the nonlinear current 
as a function of the direction of the superfluid velocity 
relative to the position of the node. This would be re- 
flected in an anisotropy of a term in the inverse penetra- 
tion depth which is linear in the magnetic field H. Early 
experimental work did not verify these predictions Q and 
it was suggested that impurity scattering 2] or nonlo- 
cal effects Q may be responsible. However, a more re- 
cent reanalysis of experiment has claimed to confirm 
the predictions An alternative proposal was given by 
Dahm and Scalapino 6] who examined the quadratic term 
in the magnetic response of the penetration depth, which 
shows a I IT dependence at low T as first discussed by 
Xu et al.Q. Dahm and Scalapino demonstrated that 
this upturn would provide a clear and unique signature 
of the nodes in the d-wave gap and that this feature could 
be measured directly via microwave intermodulation ef- 
fects. Indeed, experimental verification of this has been 
obtained^ confirming that nonlinear microwave current 
response can be used as a sensitive probe of issues as- 
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sociated with the order parameter symmetry. Thus, we 
are led to consider further cases of gap anisotropy and 
turn our attention to the two-band superconductor MgB2 
which is already under scrutiny for possible applications, 
including passive microwave filter technology 8j . MgB 2 
was discovered in 20010 and since this time an enor- 
mous scientific effort has focused on this material. On the 
basis of the evidence that is available, it is now thought 
that this material may be our best candidate for a classic 
two-band electron-phonon superconductor, with s-wave 
pairing in each channel 10] . A heightened interest in 
two-band superconductivity has led to claims of possible 
two-band effects in many other materials, both old|ll| 
and new[T^. 

Our goal is to compare in detail the differences between 
one-band and two-band s-wave superconductors in terms 
of their nonlinear response, that would be measured in 
the coefficients defined by Xu et al.Q and Dahm and 
Scalapino [13. This leads us to reconsider the one-band 
s-wave case, where we study issues of dimensionality, im- 
purities, and strong electron-phonon coupling. We find 
new effects due to strong-coupling at both high and low 
T. We then examine the situation for two-band supercon- 
ductors, starting from a case of highly decoupled bands. 
Here, we are looking for signatures of the low energy scale 
due to the smaller gap, the effect of integration of the 
bands, and the response to inter- and intraband impu- 
rities. Unusual behavior exists distinctly different from 
the one-band case and not necessarily understood as a su- 
perposition of two separate superconductors. Finally, we 
return to the case of MgB2 which was studied previously 
via a more approximate approach^|. In the current 
work, we are able to use the complete microscopic theory 
with the parameters and the electron-phonon spectral 
functions taken from band structure [Hi 13 ■ I n this way, 
we provide more detailed predictions for the nonlinear 
coefficient of MgP>2. 

In Section II, we briefly summarize the necessary the- 
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ory for calculating the gap and renormalization function 
in two-band superconductors, from which the current 
as a function of the superfluid velocity v s is then de- 
rived. In Section III, we explain our procedure for ex- 
tracting the temperature-dependent nonlinear term from 
the current and we examine the characteristic features 
for one-band superconductors in light of issues of dimen- 
sionality, impurity scattering and strong coupling. Sec- 
tion IV presents the results of two-band superconductors 
and simple formulas are given for limiting cases which 
aid in illuminating the effects of anisotropy. The case 
of MgB2 is also discussed. We form our conclusions in 
Section V. 



II. THEORY 

The superfluid current has been considered theoret- 
ically in the past by many authors for s-wavep^l Il6l 
IH El El Hi EJ Hi and for other order parameters, 
such as d-wave and f-waveQ, H HJ HI Most re- 

cently, the case of two-band superconductivity has been 
examined[ll |H H3 with good agreement obtained be- 
tween theory and experiment for the temperature depen- 
dence of the critical current 26] . 

In this work, we wish to calculate the superfluid cur- 
rent as a function of superfluid velocity v s or momentum 
q s and extract from this the nonlinear term. To do this, 
we choose to evaluate the expression for the superfluid 
current density j s that is written on the imaginary axis 
in terms of Matsubara quantities. 19, 26] This naturally 
allows for the inclusion of impurity scattering and strong 
electron-phonon coupling in a numerically efficient man- 
ner. Written in general for two-bands having a current 
j s i and j S 2, for the first and second band, respectively, 
we have: 
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where e is the electric charge, m is the electron mass, T 
is the temperature, s; = vfiIs, ni is the electron density 
and vfi is the Fermi velocity of the Z'th band (I = 1,2). 
The (•••); represents an integration for the Tth band 
which is given as J_ 1 dz/2 for a 3d band and J Q * d9/(2n) 
for a 2d band, with z = cos 9 in the 2d case. Also, in 
the expression for the current, the 3 should be changed 
to a 2 for 2d. This is done within a mean-field treatment 
and ignoring critical fluctuations near T c . Here, we have 
taken the approximation of a spherical Fermi surface in 
3d and a cylindrical one in 2d as we will see further on 
that the differences between 2d and 3d are not significant 
to more than a overall numerical factor and so providing 
more precise Fermi surface averages will not changes the 
results in a meaningful way. To evaluate this expression, 
we require the solution of the standard s-wave Eliash- 
berg equations for the renormalized gaps and frequencies 



Ai(n) = Zi(n)Ai(n) and Qi{n) = Zi(n)u> n , respectively. 
These have been generalized to two bands and must also 
include the effect of the current through q s . With further 
details given in Refs. [Toll^ . we merely state them here: 

Ai(n) = 7rTyy[Xij(m - n) - n*j{w c )Q{u) c ~ \u m \)} 



m j 



Aj (to) 



(u>j(m) — isjz) 2 + A|(m) / j 



r'Enl , Aj( '" - ) 0) 



and 



uji (n) = u n + 7rT y y Xij (to - n) 

m j 

I LUj(m) — iSjZ 



(u>j(m) — isjz) 2 + A 2 (m) / j 
u>j(n) — isjz 
fain) ~ isjzf + A 2 An) 




.(3) 



where j sums over the number of bands and the sum 



over to is from — oo to oo. Here, t 



l/(2*r+) 



the ordinary impurity scattering rate and n indexes the 
n'th Matsubara frequency u> n , with uj n = (2n — 1)ttT, 
where n = 0, ±1, ±2, •••• The /i^ are Coulomb repul- 
sions, which require a high energy cutoff u> c , taken to be 
about six to ten times the maximum phonon frequency, 
and the electron-phonon interaction enters through 



Xij (to - n) 



Qa 2 Fij (ft) 

o ft 2 + (uj n - uj m y 



-dft. 



(4) 



with ctyi^ft) the electron-phonon spectral functions and 
ft the phonon energy. Note that the dimensionality does 
not change the gap equations when there is no current. 
For finite q s , it does and we will see later the result of 
this effect. Likewise, an essential ingredient is that the 
current enters the Eliashberg equations and provides the 
bulk of the nonlinear effect for temperatures above T ~ 
0.5T C . Indeed, at T c all of the nonlinearity arises from 
the gap. 

We now proceed to the case of one-band superconduc- 
tors, to illustrate the generic features of the superfluid 
current and demonstrate how we extract the nonlinear 
term. In the section following, we will return to the two- 
band case. 



III. ONE-BAND S-WAVE SUPERCONDUCTORS 

In Fig. ^ we illustrate that these equations reproduce 
the standard results for j s versus q s for a one-band su- 
perconductor in the weak coupling BCS limit. Equa- 
tions Q)-© were solved for both the 2d and 3d cases 
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FIG. 1: The normalized current j s /jo as a function of 
Qs^f/Ao, where jo = neAo/(mtij?) and Ao is the energy gap 
at T = 0. Shown are the low temperature BCS curves for two 
dimensions (solid line) and three dimensions (dashed), given 
for a reduced temperature t = T/T c — 0.1. The inset is for 
near T c at t = 0.95. 



at t = T/Tr = 0.1 and 0.95. The T — result of past 
literature [Til [l(| is recovered in the case of 3d. One 
sees for t = 0.1 at low q s , the curve is essentially linear, 
reflecting the relationship of j s — n s ev s , with n s the su- 
perfiuid density. For strong coupling, the slope would be 
reduced by approximately 1 + A as the superfluid conden- 
sate is also reduced by this factor. Likewise the reduction 
in the slope with temperature would reflect the temper- 
ature dependence of the superfluid density. Indeed, to 
provide these curves using the Eliashberg equations, we 
used the a 2 F{u>) spectrum of Al and made the corrections 
for the 1 + A factor. Al is a classic BCS weak coupling su- 
perconductor, that agrees with BCS in every way and is 
generally used for BCS tests of the Eliashberg equations. 
The A for Al is 0.43. While at low T the curves show little 
deviation from linearity at low q s , and thus the nonlin- 
ear correction will be essentially zero (exponentially so 
with temperature in BCS theory), at T near T c , one sees 
that there is more curvature for q s — > and hence a 
larger nonlinear term is expected. However, while the 2d 
and 3d curves differ in behavior near the maximum in 
j s , one finds that the behavior at low q s is very similar. 
Indeed, the nonlincarity is a very small effect on these 
plots and hard to discern, however, it will be borne out 
in our paper that the nonlinear current does not show 
significant differences in the T-dependence between 2d 
and 3d. Nevertheless, we will still include both the 2d 
and 3d calculation in our two-band calculations as MgB2 
has a 2d er-band and a 3d 7r-band, and there is a overall 
factor of 2 /3 between the two in the nonlinear term due 
to dimensionality 

To obtain the nonlinear current as q s — > 0, the general 



expression for the current can be expanded to second 
lowest order in powers of q s leading to the general formula 
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where only first and third order terms arise. Here by 
choice jo = neAo/(mv F ) and the variable for the expan- 
sion was taken as q s v F /Ao, where v F = v F /(l + A). A 
and B are temperature-dependent coefficients which fol- 
low when solutions of the Eliashberg equations J2J an d 
(0 are substituted in the expression (|TJ for the current. 
In practice, it is complicated to expand Eqs. Q~© to 
obtain an explicit form B and so we chose to extract A 
and B numerically by solving our full set of equations 
with no approximations for j s versus q s . From this nu- 
merical data, we find the intercept and slope of j s /q s 
versus q% for q s — > from which we obtain the A and B, 
respectively. 
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FIG. 2: The nonlinear coefficient B(T) defined in Eq. © 
multiplied by (T c /Ao) 2 as a function of T/T c . The solid curve 
is for 2d BCS, while the dashed is for 3d BCS. The dot-short- 
dashed curve is for Pb in 3d, which has been evaluated within 
Eliashberg theory, and the dot-long-dashed curve shows the 
effect of impurities on Pb, where t + = T c q. The inset shows 
the temperature dependence of the linear coefficient A(T) of 
Eq. 0, which is proportional to the superfluid density. 

The results for A(T) and B(T) as a function of tem- 
perature are shown in Fig.|3 One sees, in the inset, A(T) 
which, in the one-band case, is just the superfluid density 
n s normalized to the clean BCS value at T = 0. There is 
no difference between 2d and 3d BCS. Also, shown is the 
A{T) extracted for the strong electron-phonon coupling 
superconductor Pb with no impurities and with impu- 
rity scattering of t + = T c0 . One sees that strong cou- 
pling pushes the temperature dependence of the curve 
higher, even slightly so at T = 0, and this is a well- 
documented effect |28j. With impurities, the superfluid 
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density is reduced in accordance with standard theory. 
These curves were obtained from our j s calculations and 
agree exactly with BCS and Eliashberg calculations done 
with the standard penetration depth formulas [28], con- 
firming that our numerical procedure is accurate. The 
second term in Eq. gives the nonlinear current and 
the coefficient B(T), which is a measure of this, is also 
shown for the four cases. [Note that B(T) is the same 
as the (i{T) of Ref. to within a constant of propor- 
tionality] Here one does find a difference between the 2d 
and 3d BCS curves showing that dimensionality can af- 
fect the nonlinear current. In the case of strong coupling 
one finds an increase in the nonlinear piece near T c and 
also a finite contribution at low T which is unexpected 
in the usual BCS scenario. Impurities have the effect of 
further increasing the low T contribution and reducing 
the curve near T c . 

Near T equal to T c (i = T/T c ) in BCS, it can be shown 
analytically that 



A = 2(1 - t) 



and for 3d 



c(sl a2 



(6) 



(7) 



as obtained in our previous paper, Ref. |26(. The value of 
B for 2d is increased by a factor of 3/2. These numbers 
agree with the numerical calculations in Fig. |3 where the 
3d BCS curve goes to 0.21 for 2d and 0.14 in 3d. 

There are two definitions in the literature for the 
nonlinear coefficient: one is denoted as a(T) due to 
Xu et al.J3 and the other, b(T), used by Dahm and 
Scalapino[lij, is the one that is related to the intermodu- 
lation power in microstrip resonators. Rewriting Eq. (J5J 
in the form 



, <lsV F . 



A 3 U 



(8) 



(9) 



Dahm and Scalapino define 13] 

Xu, Yip and SaulsQ keep the form of Eq. JSJ but define 
a variable q s Vp/ Aq(T), where Aq(T) is the temperature 
dependent gap equal to Ao8(t), With this they identify 
the coefficient 



a(T) = ^S 2 (t). 



(10) 



In this work, we always take S(t) to be the usual BCS 
temperature dependence of the gap function. 

In Fig. 13 we show the calculations for the a(T) co- 
efficient of Xu et al.. Here, we have made a number of 
points. First, the 2d BCS curve derived from our proce- 
dure agrees with that shown by Xu et al.0, once again 
validating our numerical work for extracting the very tiny 
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FIG. 3: The nonlinear coefficient a(T) defined in Xu, Yip, and 
Sauls0 (Eq. versus T/T c . In the top frame, we show 

curves for 2d and 3d BCS. Multiplying the 3d BCS curve by 
3/2 gives the dotted curve. The long-dashed curve illustrates, 
for the 3d case, the effect of not including the g s -dependence in 
the gap equations. In the bottom frame, the strong coupling 
Eliashberg theory result for Pb is shown, with the 3d BCS 
curve repeated for reference. The dot-short-dashed curve is 
for pure Pb and the dot-long-dashed is for t + = T c q. 



nonlinear coefficient. Second, for BCS one sees a differ- 
ence between 2d and 3d in the nonlinear coefficient. The 
2d curve goes to 1 at T c and to 2/3 for the 3d case. The 
question arises as to whether the difference between 2d 
and 3d is simply a numerical factor and so with the dot- 
ted curve, we show the 3d case scaled up by 3/2. We do 
note that there is a small difference in the temperature 
variations at an intermediate range of T, but the ma- 
jor difference between 2d and 3d is the overall numerical 
factor of 2/3. Third, one might question the necessity of 
including the effect of the current on the gap itself and 
to answer this, we show the long-dashed curve where the 
q s dependence was omitted in the Eliashberg equations 
(J2J and © . One finds that the nonlinear coefficient is re- 
duced substantially at temperatures above ~ 0.5T C and 
disappears at T c . Thus, without the q s dependence in 
the gap, the true nonlinear effects will not be obtained 
for high temperatures as the gap provides the major con- 
tribution to the nonlinearity. 

In the lower frame of Fig. [21 we examine the case of Pb 
to illustrate strong electron-phonon coupling and impu- 
rity effects. It is seen that the strong coupling increases 
the value at T c and also gives a finite value at low T. 
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The behavior at low T is surprising in light of the BCS 
result 0, but is related to the inelastic electron-phonon 
scattering which appears to increase the nonlinear coef- 
ficient at small T in a similar way to what is already 
known about the effect of impurities in BCS 13] . The 
strong coupling behavior near T c is similar to that seen 
for other quantities such as the specific heat [2^], where 
the downward BCS curvature is now turned concave up- 
ward to higher values at T c . Impurities have the effect 
of reducing the nonlinearity near T c and increasing it at 
low T. 

Once again, in BCS we can provide some analytic re- 
sults near and at T c for a(t) which provide a useful check 
on our numerical work. For three dimensions near T c : 



«(*) 



_ 7C(3) /Ap 
12^ 2 I T c 



i -t 



and, upon substituting for 5(t), 



a(t = 1) = -. 



(11) 



(12) 



Doing the same algebra for the two-dimensional case cor- 
rects these expressions by a factor of 3/2 and gives 1 
instead of 2/3 for a(t = 1). 

To characterize the strong-coupling effects seen in the 
figure for Pb, we can develop a strong-coupling correc- 
tion formula for a(t = 0) and a(t = 1). These formulas 
have been provided in the past for many quantities and 
form a useful tool for experimentalists and others to esti- 
mate the strong coupling corrections. [28| This was done 
by evaluating this quantity for ten superconductors using 
their known a 2 F(uj) spectra and their T c values. We used 
Al, V, Sn, In, Nb, V 3 Ga, Nb 3 Ge, Pb, Pb . 8 Bi .2, and 
Pbo.65Bio.35- These materials were chosen to span the 
range of typical s-wave superconductors with strong cou- 
pling parameter T c /ui\ a ranging from 0.004 to 0.2. The 
details of these materials and references for the spectra 
may be found in the review by Carbotte|28j. The param- 
eter win is defined as: 



u>l n = exp 



(13) 



By fitting to these materials, we arrived at the following 
strong coupling correction formulas for three dimensions: 



a(T c ) 



1 + 7.7 



Win 



In 



3^1n 



and 



a{T = 0) = 1.61 — 

\ Win 



(14) 



(15) 



Note that, even though ax 2 \n(b/x) is the usual form of 
the strong coupling correction, in this last equation, we 
have found no advantage in fitting with the additional 



parameter offered by the log factor. These formula should 
be seen as approximate tools to give the trend for T c j wi n 
for values restricted to the range of to 0.2. Pb has 
T c /LU\ n value of 0.128 and is intermediate to this range, 
and Al is a weak coupling superconductor with a value 
of 0.004. 
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FIG. 4: The temperature dependence of the nonlinear co- 
efficient b(T) defined by Dahm and Scalapino^| (Eq. ©) 
versus T/T c . Shown are results for 2d and 3d BCS (solid and 
dashed curves, respectively) and 3d Pb in Eliashberg theory 
(dot-short-dashed). As before these curves are for the pure 
case and one impure case for Pb is shown (dot-long-dashed 
curve) with t + = T c q. The inset shows b(T) versus T/T c , 
which is defined in Eq. 1)16^ , 

In Fig. |1 we show the coefficient used by Dahm and 
Scalapino[P3j for the same cases as previously consid- 
ered. With this coefficient one finds qualitatively similar 
curves. The 2d and 3d BCS curves go to zero rapidly 
at low temperature, but once again the strong coupling 
effects in Pb give a finite value for b(T) at low T. With 
impurities the tail at low temperature is raised signifi- 
cantly. Due to the divergence in b(T) near T c because 
of the division by three powers of the superfluid density 
which is going to zero at T c , we prefer to work with a 
new quantity b(t), which removes this divergence. Thus, 
we define 



b( t ) = b(t)(l-tf 



(16) 



and this is shown in the inset in Fig. 0] It has the ad- 
vantage of illustrating the detailed differences between 
the curves more clearly and providing finite values at T c 
which can be evaluated analytically in BCS theory. In 
this instance, we obtain 



lim b(t) 



7C(3) /Ap 
48^ 2 I T c 



0.0557 



(17) 



for three dimensions in agreement with what we obtain 
from our numerical work, shown in the Fig. 0] Once 
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again we can develop strong coupling formulas for this 
quantity and they are given as: 



b(T c ) = 0.0557 



1 - 42. 



21 

Win 



hi 



and 



b(T = 0) = b(T = 0) = 1.4 



Win 

3.8T C 



Win 



(18) 



(19) 



for three dimensions. Once again, there was no extra 
advantage to fitting b(T — 0) with the usual form that 
includes the log factor. 

This last quantity b(t) is related to the intermodula- 
tion power in microstrip resonators and hence can be 
measured directly. Having identified the features of one- 
band superconductors, we now turn to the two-band case 
where signatures of the two-band nature may occur in 
these nonlinear coefficients. 



IV. TWO-BAND S-WAVE 
SUPERCONDUCTORS 

The generalization of Eq. (JSJ) to the two-band case pro- 
ceeds as follows. The total current j s is the sum of the 
two partial currents j s i, i = 1,2 with (1,2) = (er, it) for 
the two-dimensional a- and three-dimensional 7r-band, 
respectively. For our numerical work, we do take into 
account the different dimensionality of the bands but, 
for simplicity in our analytic work below, we take them 
both to be three dimensional. A decision needs to be 
taken about the normalization of the current j s in the 
second term. Dahm and Scalapino have used jo-rr = 
n T eAo T /(m«F?r)- Here instead, we prefer to use the more 
symmetric form 



Joo 



2 



i=l 



mvFi 



(20) 



(q s /2S(t)) 2i=i v Fi/^oi: with S(t) the usual temperature 
profile of the BCS gap. Other choices could be made. 
The superfluid density n s is proportional to A for the 
combined system, specifically en s /m is given by Eq. I|21|) 
with the first two factors omitted. 
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FIG. 5: The temperature dependence of nonlinear coefficient 
B(T) (top frame) and b{T) (bottom frame) for a model two- 
band superconductor based on a Lorentzian model for the 
spectral densities afjF(w), as described in the text. Shown 
is the effect of increased interband coupling, beginning with 
a nearly decoupled curve with A12 = A21 = 0.0001 (solid) 
and progressing to more interband coupling with A12 = A21 = 
0.01 (dashed) and 0.1 (dot-dashed). The dotted curve is for 
A12 = A21 = 0.0001 and with interband impurity scattering 
of tf 2 = £21 ~ O.lTco included. 



which reduces properly to the one-band case when our 
two bands are taken to be identical, with rij = n% = n/2, 
where n is the total electron density per unit volume. For 
the combined system, Eqs. (JHJ and (|10|1 still hold with A 
and B modified as follows: 
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2 A 0l 
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and 
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Fi 
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(22) 



With these definitions Eq. (JSJ also holds with j o 
replacing jo and the Xu, Yip, and Sauls variable, 
q s v F /A (T), of the one-band case is replaced by 



In Fig. we show both the B{T) and the b(T) for a 
model which uses truncated Lorentzians for the a 2 Fij(ui) 
spectra. This same model was used in our previous 
work^ an d so we refer the reader to that paper for 
details. Also, in Ref. 01 mav be found the curves 
for the Aj(T), the penetration depth, and other quan- 
tities for the same parameters used here. The essen- 
tial parameters of this model are An = 1, A22 = 0.5 
and the interband electron-phonon coupling is varied 
from A12 = A21 = 0.0001 (nearly decoupled case) to 
0.1 (more integrated case). In addition, the = 0, 
Vpi = vp2 and n\ — 712. In the nearly decoupled case of 
A12 = A 2 i = 0.0001, it can be seen that the solid curve 
looks like a superposition of two separate superconduc- 
tors, one with a T c which is about 0.33 of the bulk T c . 
The lower temperature part of this curve is primarily due 
to the 7r-band (or band 2) which is three dimensional, and 
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indeed, when examined in detail, it has the characteris- 
tic behavior of the 3d example studied in the one-band 
case. The part of the curve at higher temperatures above 
about t ~ 0.35 is due to the cr-band (or band 1) which is 
taken to be 2d and indeed, in the case of B(T) it shows a 
dependence approaching T c that expected for 2d strong- 
coupling with some interband anisotropy effects. The rel- 
ative scale of the two sections of the curve is set by the 
value of the gap anisotropy u = A 2 / A i and the ratio 
(1 + An +Ai 2 )/(1 + A 2 2 +A 2 i). The overall scale on the y- 
axis for b(T) differs from that of Fig. 0] due to our choice 
°f Joo for the normalization in the nonlinear term. In- 
deed, for nearly decoupled bands (solid curve), the value 
of the nonlinear coefficient b(T) is small at reduced tem- 
perature t = 0.4 just above the sharp peak due to band 2. 
Specifically, it is of order 0.5. If it had been referred to jo a 
instead of joo, it would be smaller still by a factor of 1.7 
and comparable to the single band 2d BCS result at the 
same reduced temperature (Fig. 0] bottom frame, solid 
curve). However, as the non-diagonal electron-phonon 
couplings A 12 and A 2 i are increased and a better inte- 
gration of two bands proceeds, b(T) at t = 0.4 can in- 
crease by an order of magnitude as, for example, in the 
dashed curve. The actual scale in this region is set by 
the details of the electron-phonon coupling (see later the 
specific case of MgB 2 ). 

With more integration between the bands, one finds 
that the sharp peak at lower T is reduced and rounded 
with a tail reaching to T c . When A12 = A 2 i = 0.1, the 
feature characteristic of the 7r-band T c is almost gone in 
B(T) and absent entirely in b(T), even for modest inter- 
band coupling. The same conclusion holds for the effects 
of interband scattering (shown in Fig. [3] for a value of 
<+ = t+ = 0.IT c0 for the nearly decoupled case) which 
also integrates the bands and eliminates the lower en- 
ergy scale. However, while the structure at the lower T c 
is now reduced to the point of giving a monotonic curve 
for &(T), there still remains a large nonlinear contribu- 
tion well above that for the one-band s-wave case, which 
marks the presence of the second band. 

We can have further insight into these results and check 
our work by developing some simple analytic results in 
renormalizcd BCS theory (RBCS). For a summary of the 
approximations of RBCS and a comparison with full nu- 
merical solution for various properties including j s , we 
refer the reader to our previous work[l(l|2(l. For sim- 
plicity, we take both bands to be three dimensional in the 
following. Near T = T c , Eqs. © and are modified for 
each band to: 



Ai = 2(1 - t) — 

Xi 



and 



Bi 



7 C(3) 



A 2 

^Oi 



6(ttT c ) 2 X \ 



(23) 



(24) 



where the functions \i ancl Xi have been derived in 
Ref. [2^| and v F \x'i is independent of v Fi , where v* Fi = 



VFi/{l + The x's depend on the microscopic 

parameters of the theory. In RBCS, they are Ay, fi*p 
Vpi, and rii, from which T c and Aoi follow. While 
the expressions obtained for the \i an d Xi are lengthy, 
and hence we do not repeat them here, they are ex- 
plicit algebraic forms. It is useful in this work to con- 
sider several simplifying limits. For decoupled bands 
A12 = A 2 i = /i* 2 = M21 = an d A 2 = B 2 = 0. As the 
band 2 does not contribute near T c , A\ and B\ take on 
the form of the single band case (Eqs. © and Q). An- 
other limiting case is the separable anisotropy model. j2|J 
In this model, there are only two gap values with a ratio 
of u = (1 — a) /(I + a), with a an anisotropy parameter 
often assumed small. In this case, An = A(l + a) 2 /2, 
A 22 = A(l - a) 2 /2 and A i2 = A21 = A(l - a 2 )/2, where 
A = A/(I + A) and Ay = Ay/ (I + J2i ^0- As a resu lt 



1 

Xi 
and 



(1 + a) 2 (l - 5a 2 ), — = (1 - a) 2 (l - 5a 2 ), (25) 

X2 



Xi Xi 



(26) 



In this model, taking in addition that vfi — v F 2 and 
n± = n 2 = n/2 leads to the one-band case and this can 
be used as a check of our algebra. 

To see the consequences of this algebra for our nonlin- 
ear coefficient b(T), we begin with the decoupled band 
case near T = T c for which ^4 2 = B 2 = and A\ and B\ 
reduce to their single band value. In this limit of t — > 1, 



m = I 



i 



I 

VF2 



(27) 



where u is the gap anisotropy parameter u = A 02 /A i, 
and Aoi is the gap at T = 0. This expression shows ex- 
plicitly the corrections introduced by the two-band na- 
ture of the system over the pure one-band case. Note 
that b(T) is always increased by the presence of the cor- 
rection term. In l|27|) . u can never be taken to be one 
since we have assumed band 2 is weaker than band 1. 
Before leaving the decoupled case, it is worth noting that 
b(T) will show a change at the band 2 critical temper- 
ature T c2 . For T below T C 2, A 2 and B 2 will be finite 
while above this temperature they are both zero. When 
the coupling Ai 2 and A 2 i is switched on but still small, 
we expect that these quantities will acquire small tails 
and that they vanish only at T c . This is the hallmark of 
nearly decoupled bands. 

For the anisotropic a 2 model near T c 



b(t) 



_ 7C(3) /Ag"\ 2 I + 8a 



6ir 2 \T C J [2(1 - t)} 



(28) 



where the average gap Aq" is related to T c by 2Aq"/T c — 
3.54[1 — 3a 2 /2]. For a 2 = this expression reduces prop- 
erly to the one-band limit. Therefore, it is seen that 
anisotropy increases b(T) for T near T c . 
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Another interesting limiting case is to assume both 
bands are the same, i.e. isotropic gap case, but that 
the Fermi velocities differ in the two bands. Near T c , we 
obtain 



1 



1 (VFI 



VF2) 



7C(3) (A,/ 

6tt 2 \T C J [2(1 -t)] s 8 {v Fl v F2 f 



(29) 

In this case, the Fermi velocity anisotropy changes the 
nonlinear coefficient, but when Vfi = v F 2 the expression 
reduces properly to the one-band result. We find that 
the Fermi velocity anisotropy increases b(T) near T c , a 
result that is seen in one of our calculations for MgB2 
shown in Fig. |SJ 
0.6 




0.0 0.2 0.4 0.6 0.8 1.0 
T/T 




FIG. 6: The nonlinear coefficients B(T) (top frame) and b(T) 
(bottom frame) versus T/T c for the case of MgE>2. The solid 
curve is for the pure case (no impurity scattering), the dashed 
for impurities in the a band with tf x — T c o and dot-dashed for 
tj 2 = Tco- In the pure case, increasing the ratio of vf2/vfi = 
3 gives the dotted curve. The inset in the top frame gives 
A(T) vs T/T c and the bottom frame inset shows b(T). 

With Fig. [HI we now turn to the specific case of MgB2, 
where we have used the parameters and a^-F^w) given 
by band structure calculations, and as a result, there 
are, in principle, no free parameters other than varying 
the impurity scattering rate. The basic parameters are 
A (j (j = 1.017, = 0.448, X a7T = 0.213, K a = 0.155, 
M<7<T = 0.210, = 0.172, ^ = 0.095, ^ a = 0.069, 
with aT c — 39.5 K and a gap anisotropy of u = 0.37. The 
ratio of the two density of states is A^OyA^O) = 1.37 
and of the Fermi velocities is v F ^jv F(J — 1-2. We have 
found excellent agreement between theory and experi- 
ment for these parameters, as have other authors [Tol l26|| . 



As we have found in our previous work, MgB2 is quite 
integrated between the bands. It is also an intermedi- 
ate strong coupler with T c /uj\ n — 0.05 and thus there 
is competition between the strong coupling effects and 
the anisotropy TOJ . In Fig. 6, the solid curve gives the 
prediction for MgB2 for B{T) and b(T). A strong non- 
monotonic feature around the lower band energy scale is 
observed in B(T), but the b(T) is monotonically increas- 
ing with temperature. To see the second band effects in 
b(T), it is better to plot b(T) (the inset) which accentu- 
ates the subtle variations found at the lower energy scale 
associated with the tt band. Also shown in the inset for 
the upper frame is A(T), which gives the temperature 
dependence of the superfluid density. The solid curve 
agrees with our previous calculation by other means [T(l|. 
The variation in A(T) appears to be sufficient to remove 
the bump in B(T) when divided by three factors of A(T) 
to obtain the definition of b(T). Also, shown are the ef- 
fects of intraband scattering with = T c0 for the dashed 
curve and i J 2 = T c o for the dot-dashed. Scattering in the 
•rr-band reduces its contribution and provides an impurity 
tail at low T, as found for the one-band case. However, 
scattering in the cr-band, while lowering B(T) near T c as 
expected, does not appear to add weight at low T. This 
is because the parameters for MgB 2 heavily weight the w- 
band and the cr-band is a small component. Thus, upon 
comparison between cr- and 7r-band scattering, b(T) could 
be lowered at t — 0.5, for example, by putting impurities 
in the 7r-band, but it would be raised if the impurity scat- 
tering is in the cr-band. The dotted curve in the figure 
is for pure MgB 2 but where we have taken v F2 /v F1 = 3 
to mimic a case where transport may happen along the 
c-axis. In this instance, the bump in B(T) remains, but 
is gone in b(T). We see that b(T) is large due to the 
higher power of the Fermi velocity ratio that enters the 
calculation, and, as a result, the nonlinearity is greatly 
increased. As b(T) is a relevant quantity for microwave 
filter design, this study provides some insight into which 
factors may be used to optimize the material and reduce 
the nonlinear effects. 



V. CONCLUSIONS 

Study of nonlinear current response is important for 
device applications and for providing fundamental sig- 
natures of order parameter symmetry, such as have been 
examined in the case of d-wave superconductivity. In this 
paper, we have considered the case of two-band supercon- 
ductors. In so doing, we also reexamined the one-band 
case and discovered that there can exist extra nonlin- 
earity at both low and high temperatures due to strong 
electron-phonon coupling, for which we have provided 
strong-coupling correction formulas, whereas the excess 
nonlinearity induced by impurity scattering occurs pri- 
marily at low temperatures. At T c , impurities will give 
an enhanced or decreased contribution depending on the 
particular nonlinear coefficient discussed. In this paper, 
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we have examined two nonlinear coefficients denned in 
the literature, one due to Xu et al.0 and one defined by 
Dahm and Scalapino, with an emphasis on the latter as 
it is related to the intermodulation power in microstrip 
resonators jl3|. 

We have also studied issues associated with dimen- 
sionality motivated by the two-band superconductor 
MgB2 , which has a two-dimensional er-band and a three- 
dimensional 7r-band. Within our one-band calculation, 
aside from an overall factor of 2/3, we find little differ- 
ence in the temperature variation of the nonlinear coef- 
ficient in mean-held between 2d and 3d. This is further 
reduced by strong coupling effects. 

For two-band superconductors, we show that for nearly 
decoupled bands a strong signature of the small gap 
7r-band will appear in the nonlinear coefficients, but 
with increased interband coupling or interband scatter- 
ing, such a signature will rapidly disappear. Likewise, 
intraband impurities in the 7r-band will wash out the tem- 
perature variation of the 7r-band, whereas the intraband 
impurities in the er-band largely effect the nonlinearity 
at higher temperatures above the energy scale of the w- 
band, for the parameters typical to MgB2. 

We provide a prediction for the nonlinear coefficient in 
MgB2 using the parameters set by band structure calcu- 
lations. As the bands in MgB2 are quite integrated, we 
find that the nonlinear coefficient b(T) is monotonically 
increasing in contrast to a previous prediction, which was 



based on a number of approximations, |f3j and we find 
that the increased nonlinearity due to the 7r-band is best 
reduced at t = 0.5 by adding impurities to the 7r-band. 
Should the supercurrent sample the c-axis direction, a 
larger anisotropy in the Fermi velocity ratio between the 
bands would result and this effect is found to increase 
the nonlinearity. Finally, several simple formulas have 
been provided for near T c which aid in the understand- 
ing of the range of behavior observed in the numerical 
calculations. We await experimental verification of our 
predictions. 
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